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Abstract: In this paper, we study a Lotka-Volterra model which contains two prey and one predator with the Beddington- 
DeAngelis functional responses. First, we establish a set of sufficient conditions for existence of positive periodic so¬ 
lutions. Second, we investigate global asymptotic stability of boundary periodic solutions. Finally, we present some 
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1. INTRODUCTION 


predator consumes prey with functional responses: 


The dynamical relationship between predators and 
prey has been studied by several authors for a long time. 
In those researches, various forms of functional responses 
have been used. Here a functional response means the 
average number of prey killed per individual predator per 
unit of time. Some biologists have argued that in many 
situation, especially when predators have to search for 
food, functional responses should depend on both prey’s 
and predator’s densities, see [1,7,12,13] and references 
therein. 

Let us consider a population of three species, say a 
Lotka-Volterra model, with the following properties: 

(i) one species is a predator of two competitive other 
species. 

(ii) the predator consumes prey with the functional re¬ 
sponse given by Beddington [2] and DeAngelis et al. [ 6 ] 

There are many models having the property (i) or (ii) with 
diffusion in a constant environment [3-5], [11], [15-17], 
However, natural environments are usually periodic in 
time due to the periodicity of seasons. Therefore, the pa¬ 
rameters in these models should be periodic in time. This 
paper devotes to studying such a Lotka-Volterra model 
which is performed by a nonlinear system of differential 
equations: 
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( 1 . 1 ) 


Here 27 (t) represents the population density of species 
Xi at time t{i 1), X :i is a predator species and Xi,X 2 
are competitive prey species. At time t, a^t) is the in¬ 
trinsic growth rate of Xj (i = 1 , 2 ) and 03(f) is the death 
rate of X 3 ; bij(t) measures the amount of competition 
between X\ and X 2 (i j,i,j ^ 2), and ba(t) (i ^ 2) 
measures the inhibiting effect of environment on Xi. The 


_ Cl (t)xix 3 _ _ c 2 (t)x 2 x 3 __ 

a(t) + /3(t)xi + 7 (t)x 3 an a(t) + @(t)x 2 + 7 (f)x 3 ’ 

and contributes to its growth with amounts: 

_ d\{t)xi _ _ d 2 {t)x 2 _ 

a{t) + fd(t)x\ + j(t)x 3 a(t) + /3(t)x 2 + 'y{t)x 3 

Furthermore, we assume that the parameters a,; (f), (f), 
Cj(f), di(t), a(t), /3(f), 7(f)(l < i,j ^ 3) are w-periodic 
and continuous in t and bounded below by some positive 
constants. 

In the next section, we present our main results. First, 
we use the continuation theorem in coincidence degree 
theory to show existence of positive periodic solutions 
of (1.1). Second, by using Lyapunov functions we ver¬ 
ify global asymptotic stability of boundary periodic solu¬ 
tions. Finally, we give numerical examples. 

2. MAIN RESULTS 

For biological reasons we only consider (1.1) with 
nonnegative initial values, i.e. a 3 i( 0 ),X 2 ( 0 ),x 3 ( 0 ) ^ 0 . 
Let g(t) be a function, for a brevity, instead of writing 
g(t) we write g. If g is a bounded continuous function on 
R, we denote 

9 U = sup tgR g(t), g l = inf( g R g(t), 

and f/=J- j 0 g(t)dt, if g is a periodic function with pe¬ 
riod LO. 

Definition 1: A nonnegative solution x*(t) of (1.1) is 
called a global asymptotic stable solution if it attracts any 
other solution x(t) of ( 1 . 1 ) in the sense that 

lilllt-s-oo Ei =1 Mf) - x* (f) | = 0. 

2.1. Existence of positive periodic solutions 

In this subsection, we shall study existence of periodic 
solutions of (1.1). It is not difficult to verify global ex¬ 
istence and uniqueness of nonnegative solutions of ( 1 . 1 ). 
To show the existence of a positive periodic solution, we 
shall use the continuation theorem in coincidence degree 
theory which has been used for some mathematical mod¬ 
els of Lotka-Volterra type [10,16] and references therein. 












The following are some concepts and results taking from We now put 

[ 8 ]. 


Let X and ¥ be two Banach spaces. A linear mapping 
L: V(L) C X —>• Y is called Fredholm if it satisfies two 
conditions: 


(i) Im L is closed and has finite codimension; 

(ii) Ker L has finite dimension. 


The index of L is the integer dim Ker L — codim Im L. 
If L is Fredholm of index zero, there exist continuous 
projections P: X —> X and Q: Y —> Y such that ImP = 
Ker L, Im L = Ker Q = Im(7— Q), and an isomorphism 
J: Im Q —» Ker L. It follows that 

Lp = L\v(L)nKerP ■ (/ —P)X—>ImL 
is invertible. We denote the inverse of that map by K p . 
Let fl be an open bounded subset of X. A continuous 
mapping N: X — > Y is said to be L-compact on Cl if the 
following two conditions take place: 
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= In 


f a 1 -6 12 e ff21 -(^) 1 
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j-a 2 - b 2 ie Hl1 - (f) ^ 
l Ln y 


■ d\e Hl1 + foe 1121 — 0 , 3 a 1 1 
t 037 * /’ 

H31 = L31 + 2S3W, 

L 32 = ln[(di ^ d 3 /3 u )e Hl2 

+ (d 2 - a 3 /3 u )e H22 - 2 a 3 a u ] - ln( 2 a 3 7 M ), 
H 32 = L 32 — 2a 3 u>. 


(i) the mapping QN: Cl —> Y is continuous and 
bounded; 

(ii) K p (I — Q)N: id —>■ X is compact, i.e. it is continu¬ 
ous and K p (I — Q)N(Cl) is relatively compact. 


To introduce the definition of the degree of N in Q, for 
simplicity we assume that X = 1 N . Suppose furthermore 
that N is smooth on Cl. Let p Of! be a regular value 
of N, i.e. the equation N(x) — p on Cl has only a finite 
number of solutions x\,... ,x n £ fl with nonsingular 
DN(xi) for each i = 1,..., n where DN(xi ) is the Ja¬ 
cobi matrix of N at 27 . Then the degree deg(_ZV, Q. p) of 
N in fl at p is defined by the formula 


n 

deg (N,Cl,p) = y^sgnjdet DN(xj)}. 

i= 1 


Lemma 2 (Continuation theorem [ 8 ]) Let L be a Fred¬ 
holm mapping of index 0. Assume that N: Cl —» Y is 
L-compact on Cl and satisfies conditions: 


The convention here is that In 2 ; = —00 if x 0. Under 
the conditions 


bllb 2 2 7 ^ &12&21) 

01 - b 12 e H2 ' - (g) > 0 , 

< a 2 - b 2 ie ffl1 - (a) > 0 , 
d\e Hl 1 + d 2 e H21 — a 3 a z > 0, 
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( 2 . 1 ) 

we shall verify existence of an w-periodic solution of 

(1.1). 

Theorem 3: Let (2.1) be satisfied. Then (1.1) has 
at least one positive w-periodic solution. 

Proof: By putting 27 (f) = e Ui< d\i ^ 1), (1.1) be¬ 
comes 
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Let 


(a) for each A £ (0,1) every solution of Lx = A Nx is 
such that x ^ Of l, 

(b) QNx ^ 0 for each x £ OfI Cl KerL, and 
deg {JQN, f l D KerL, 0} ^ 0. 


X = Y 

= {it = (tti,u 2 , U 3 ) t £ C , 1 (R, R 3 ) such that 
Ui(s) = Ui(s + u>) for s £ R and i ^ 1}, 

with norm 


3 

IMI = max |itj(s)|, 
o,w] 


Then the operator equation Lx = Nx has at least one 
solution in 2 ?(L)(~I Cl. 


it € X. 



Then both X and Y are Banach spaces. Let 


interval [0, w], we obtain 


u\ Ni (s) 

u 2 (s) = N 2 {s) 

« 3 J L^3(s). 

Tai - 6ne“ 1 - b 12 e U2 - 0 + ^lf+ 7e « 
a 2 - b 2 ie ui - b 22 e U2 - Q+/3 %+ 7e , 
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Hence, 


pUJ 

Kei L = R 3 ,ImL = {u £ Y| / Ui(s)ds = 0, * ^ 1}, 

Jo 

and dimKerL = 3 = codim ImL. Then, it is easy to 
obtain the following conclusions. 

1. L is a Fredholm mapping of index zero, since Im L is 
closed in Y. 

2. P and Q are continuous projections such that Im P = 
KerL,ImL = KerQ = Im(/ — Q). 

3. The generalized inverse (to L) I\p: Im L —s- V(L) n 
Ker P exists and is given by 

ui] \fq u i ( s )ds - b So So M i ( s)dsdv 

K P u 2 (u)= f^u 2 (s)ds - ± S^S^M^dsdv . 
U 3 J L/o u 3 (s)ds - i So So u 3{s)dsdv_ 


4. QN and Kp(I — Q)N are continuous. 

5. TV is L-compact on 0 with any open bounded set fl C 

X. 

Let us now find an appropriate open, bounded subset 
0 for application of the continuation theorem. Obviously, 
from (2.1), —00 < La ^ Ln < 00 (i ^ 1). Correspond¬ 
ing to the equation Lu = XNu, A € (0,1), we have 

u i = A ai - 6 ne Ul - bi 2 e “ 2 - 1 —— , 

L a + pe Ul + je U3 J 

coe U3 

14 = A a 2 - b 2 ie Ml - b 22 e U2 -- , 

2 L a + Pe u2 + 7e“3 J 


u 3 = X — 03 


a + (3e Ul + 7 e “ 3 
d 2 e U2 

a + f3e U2 + 7 e“ 3 . 


Suppose that ( Ui,u 2 ,u 3 ) t £ X is an arbitrary solution 
of (2.3). Integrating both the hand sides of (2.3) over the 


pOJ 
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f‘Uj 

a 2 uj = / b 2 ie Ul + b 22 e U2 + 
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- 7 :- at, 
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(2.4) 


a + /Je^ 1 + 7 e tl3 a + /3e“ 2 + 7 e M3 . 


Combining the first equations of (2.3) and (2.4), we ob¬ 
serve that 
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+ / b 12 e U2 dt+ / - -2 - \dt 

Jo Jo a + (3e Ul + pe U3 J 

< 2 ai u>. 

Similarly, we have [J \u' 2 {t)\dt < 2 a 2 co, and 


pi0 pUJ 

i 3 (t)\dt < A / a 3 df + / — 
L do Jo a 

J d 2 e U2 i , 

1 --- dt < 2a,3U. 

Jo a. + /3e“ 2 + ~/e U3 J 


0 a + f3e Ul + je U3 


Since u £ X, there exist £i,rji £ [0, w] (* ^ 1) such that 

uJCi) = min Ui(t), Ui(rji) = max Ui(t). (2.5) 
te[o,u] te[o,w] 

From the first equation of (2.4) and (2.5), we obtain 

pOJ pUJ 

a iw^ / b 11 e Ul ^ l) dt+ / b 12 e U2 ^ 2) dt 

Jo Jo 

= 6nwe“ l(il) +b 12 coe U2 ^ 2 \ 

which implies that ui(^i) < Lu. Hence, for all t ^ 0 

pUJ 

u\(t) < ui(^i) + / \u'-y(t)\dt < Lu+2a,\uj = Hu. 

Jo 

Similarly, we have u 2 (t) < H 2 1 for all t ^ 0. 

On the other hand, from the first equation of (2.4) and 
(2.5), 


< [ b lie u ^dt+ f b l2 e U2(r,2) dt+ [ dt 

Jo Jo Jo 7W 

= Snwe" 1 ^ 1 ) + b 12 ue U2 ^ + (—)w 

7 

< S n we Ml(?7l) + 6i 2 we ff21 + (—)w. 
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Hence, 

pUJ 

u \( t ) ^ 1x1(771) — / ^ -^12, vt > 0. 

Jo 



Similarly, t* 2 ( 6 ) p Il-n for all 6 ^ 0 . Therefore, by 
putting Bi = max{|Ppi|, \Ha\}* we conclude that 


max \‘Ui(t)\ P Bi, i = 1,2. 

*e[o,w] 

Let us give estimates for 1*3(6). It follows from the 
third equation of ( 2 . 4 ) and ( 2 . 5 ) that 

d\(t)e Hl1 d. 2 (t)e H21 


c 13 UJ P 
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Hence, ^3(^3) P P31 and 1*3(773) ^ P32. We then ob¬ 
serve that 

7*3(6) < 7* 3 ($ 3 ) + fo l u 3(f)l rfi < H 31 

and 


pCO 

U3(t) > ** 3 ( 773 ) - / KWMf 

Jo 


> H- 


32- 


Therefore, by putting B3 = max{|if3i|, IIT32I}, we get 


max 17 * 3 ! p H 3 . 

£E[ 0 ,cd] 


By carrying out similar arguments as above, we observe 
that any solution it* of the equation fT M (it) = 0 £ R 3 
with 71 € [0,1] satisfies the estimate 

Li 2 < it* P P;i, i > 1 - ( 2 . 6 ) 

Thus, 0 ^ iT M (<9flnKerP) for 7 * G [0,1], Consequently, 
by taking 71 = 1 , we conclude that 0 ^ Q6V(<9GnKer L). 
Note that the isomorphism J can be the identity mapping 
I, since ImP = Ker P. By the invariance property of 
homotopy, we obtain that 


deg (JQN, Q (~l KerP, 0) = deg(QN, fl n KerP, 0) 
= deg(QN, fl n R 3 ,0) = deg(G, Q n R 3 ,0) 

= sgn det A 


= sgn-j (611622 - 612621 ) 

where 


d[f{ui,u 3 )+g(u 2 ,u 3 )\ 

du 2 


A = 


4ue“> 
-62 re Ul 

df(u 1 ,u 3 ) 
du 3 


-6i 2 e“ 2 0 

—6 22 e“ 2 0 

dg(u 2 ,u 3 ) df(ui,u 3 ) , dg{u 2 ,u 3 ) 

du% du$ ' dus 


Since both functions f(u\, 1 * 3 ) and <7(1*2, 7 * 3 ) increase in 
1 * 3 , d ^Qu™ 3 ' > + d9 ^du 3 3 ' > > H ence - by using the first 
condition in (2.1), we conclude that 

deg (JQN, Q n KerP, 0) ^ 0 . 

By now we have proved that fl verifies all require¬ 
ments of Lemma 2 . Therefore, the equation Lu = Nu 
has at least one solution in Z?(P) (~l Cl, i.e. (2.2) has 
at least one w-periodic solution it* in X>(P) 0 Cl. Set 
x* = e Ui (i ^ 1), then x* is an w-periodic solution of 
( 1 . 1 ) with strictly positive components. It completes the 
proof. ■ 


By the above estimates, for any solution u G X of 
( 2 . 4 ) we have ||it|| SC X^*=i Bi- Clearly, Bi[i ^ 1 ) are 
independent of A. Take B = 5^t=i where B 4 is taken 
sufficiently large such that B4 p YC,i= 1 1 \Lij\- Let 

fi = {« € X | ||t*|| < B}, then O satisfies the condition 
(a) of Lemma 2 . 

Let us verify that the condition (b) of Lemma 2 is also 
satisfied. Consider the homotopy 


2.2. Global asymptotic stability of boundary periodic 
solutions 

In this subsection, we shall establish a sufficient crite¬ 
ria for global asymptotic stability of boundary w-periodic 
solutions of ( 1 . 1 ). Consider the boundary dynamics of 
(1.1) where X3 is absent, i.e. X3 (t) = 0 for every 6^0. 
We then consider the periodic competitive model of two 
prey Xi,X 2 \ 


H^u) = 7 iQN(u ) + (1 - n)G{u), 7* G [ 0 , 1 ] 

where G: R 3 —>• R 3 , 


G(u) 


01 - 6 n e Ml - 6 12 e “ 2 
d 2 - 6 2 ie“ x - b 2 2e U2 
.-a 3 + /(«!, 7*3) + 3(7*2,7*3). 


with /(1*1,7*3) = 

J_ rM d, 2 e u3 dt 
co Jo a+/3e u 2 +-ye u 3 


1 rv d^e u 3 dt 
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and 3 ( 1 * 2 , 773 ) = 


. It is easy to see that 


H^u) 

= 0,2- hie " 1 - 6 22 e “ 2 - L / 0 “ a+ ^ c e P e 2 +* 

-a 3 + /(i*i, 1*3) + g(u 2 , 7*3) 


= xi [ai(t) - 6 n(t)ati - 612 ( 6 ) 0 : 2 ], 
x ' 2 = x 2 [a 2 (t) - 621 ( 6 ) 0:1 - 622 ( 6 ) 0 : 2 ]. 


( 2 . 7 ) 


Denote by X,(6) the unique positive ^'-periodic solution 
of the logistic equation: 


Then 26 , (6) = 


X' = X [a*(6) - b u {t)X] 

e fo a % (s)ds — 1 


/ t t+ “ Li(s)e“ 6»* a < T > eiT ds ' 


Due to [ 9 ], if 


at > bijXj (i ± j, i,j = 1,2), (2.8) 


then ( 2 . 7 ) has a positive w-periodic solution (o:i,o,’ 2 ). 
Furthermore, if 


A 12 < 0 


( 2 . 9 ) 



then (xi,X 2 ) is globally asymptotically stable, where 

dij(t) = bij(t)xj(t) (i ^ j, i,j = 1 , 2 ) and 

4 fj.\ f V" Q'ji) ■ / -1 

A 12 (t) = max |-—- ajj,i ± j | . 

Our result is as follows. 

Theorem 4: If (2.8) and (2.9) hold then x = 
(xi,X2,0) is a w-periodic boundary solution of (1.1). 
Furthermore, 

(i) If bij < bjj (1 < i ^ j < 2) and ci + C 2 + d\ + d 2 < 

pa 3 then x is globally asymptotically stable. 

(ii) If ci+C 2 + di+d ,2 < (3a 3 then x attracts any solution 
x of ( 1 . 1 ) which satisfies the condition 

[x\ (f) -xi(t)][x 2 (t) -x 2 (f)] >0, Vf Ss 0. 

(hi) If (h + d‘2 < fda 3 then x attracts any solution x of 
( 1 . 1 ) which satisfies the condition 


x»(f) ^ 27 (f), Vf ^ 0 ,i = 1,2. 

Proof: The first statement is obvious. To prove (i), 
let x be any other solution of (1.1). Consider a Lyapunov 
function defined by V(t) = JL=i I lnXj — lnxj|+x 3 ,f ^ 
0. Calculating of the right derivative D + V[t) of V(t) 
along the solutions of ( 1 . 1 ) gives 


D + V(t) 


= X] s g n ( x t - *»)(— - -r 1 ) + 4 

Xi Xi 

i—l 

2 

= ^ { sgn(xj - x t ) (a* - buXi - 6 , 




CiXiX 3 
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{- a 3 + Y 


) {ai baXi bijXj) ^ 
diXi 
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)X3 


= y^ {[-bu\xj - Xi\ - bij(xj - Xj)sgn(xi - Xj)]| 
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x 3 \ 


~ a 3 + E { 
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diXi 


a + /3xj + 7 x 3 

_ CjXj sgn(xj - Xj) ' 
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. ( 2 . 10 ) 


Then 


D + V(t) V - b,j) xj - Xj| 

(ci + c 2 + di + V 2 - /?a 3 )x 3 


( 2 . 11 ) 


By assumptions in (i) and the periodicity of parameters, 
there exist p± > 0 such that 

f Cl + c 2 + Vi + d 2 - P a 3 

max 1 - o - 



Thus, by integrating both the hand sides of (2.11) from 0 
to t, we observe that 

V{t) + pi f* EL \ x i - X M S < ^(0) < co 
for every t ^ 0. Hence, E^i I x % ~ %i\ € ^ 1 ([0;°°))- 
On the other hand, by the periodicity, xy and x, (1 7 1) 
have bounded derivatives on [0, 00 ). As a consequence, 
Ei=i I x i — Xi | is uniformly continuous on [ 0 , 00 ). There¬ 
fore, by using the Barbalat lemma [16], we conclude that 
lim X)i=t \ x i ~ x i\ = 

t—>00 

i.e. x is globally asymptotically stable. 

Similarly, we obtain the conclusions in (ii) and (iii) by 
using the following inequalities, respectively. 


D + V(t) 


= + bji)\xi - Xi\] + x 3 
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Ct3 


i=l 
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CiXi sgn(xi - Xi) y 

l a 

+ f3 Xi + 7 x 3 
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V ^ | {pa V bji ) | Xj Xi\\ 

, (ci + c 2 + di + d 2 - /3 a 3 )x 3 

P ’ 


and 


D+V{t) 


= '^2[-{bu + bji)\xi - Xi\] + x 3 
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diXi 


- a 3 


CiXi 


“ ^ a + fai + 7X3 a + fixi + 7X3 . 


^ ^ v [ (bn bji)\Xi XiW ~\~ 


(di + d 2 - pa 3 )x 3 


We complete the proof. ■ 

2.3. Numerical examples 

In this subsection, we exhibit some numerical ex¬ 
amples which show the convergence of positive so¬ 
lutions of (1.1) to periodic solutions of (1.1). Set 
ai = 3 + sin( 8 f);d 2 = 5.5 — 0.2 cos( 8 f); a 3 = 
0.4 — 0.3cos(8f); 6 n = 2 + cos( 8 f );&22 = 5 + 
0.4sin(8f); 612 = 0.04 — 0.02 sin( 8 f); 621 = 0.15 — 
0.1 cos( 8 f); Ci = 0.5 — 0.4sin(8f); C 2 = 0.4 — 
0.3sin(8f); a = 0.03 — 0.02 cos( 8 f); /? = 0.3 + 
0 . 2 cos( 8 f );7 = 2 — sin( 8 f);di = 3 + 2 sin( 8 f);d 2 = 
3 — 2sin(8f); and an initial value (xi(0), X 2 ( 0 ), x 3 (0)) = 
(0.5,0.7,1). Figure 1 shows the behavior of the solution 
of (1.1). It is seen that the solution converges to a positive 
periodic solution of ( 1 . 1 ). 

We now set a 3 = 4 — 0.3cos(8f) and (3 = 3 + 
0.2 cos( 8 f) and retain other parameters as above. Figure 
2 gives the behavior of the positive solution of (1.1). It 
converges to the boundary periodic solution of ( 1 . 1 ). 


< -fti- 





Fig. 1 A solution of (1.1) which converges to a positive 
periodic solution of (1.1) 



Fig. 2 A positive solution of (1.1) which converges to the 
boundary periodic solution of (1.1) 
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